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ABSTRACT 

We report the detection of seven low mass companions to intermediate-mass stars (SpT B/A/F; 
M~1.5-4.5Mq) in the Scorpius-Centaurus Association using nonredundant aperture masking inter¬ 
ferometry. Our newly detected objects have contrasts AL'«4-6, corresponding to masses as low as 
~20 Mj up and mass ratios of g~0.01-0.08, depending on the assumed age of the target stars. With 
projected separations p~10-30AU, our aperture masking detections sample an orbital region previ¬ 
ously unprobed by conventional adaptive optics imaging of intermediate mass Scorpius-Centaurus 
stars covering much larger orbital radii (~30-3000 AU). At such orbital separations, these objects 
resemble higher mass versions of the directly imaged planetary mass companions to the 10-30 Myr, 
intermediate-mass stars HR 8799, (3 Pictoris, and HD 95086. These newly discovered companions span 
the brown dwarf desert, and their masses and orbital radii provide a new constraint on models of the 
formation of low-mass stellar and substcllar companions to intermediate-mass stars. 

Subject headings: instrumentation: adaptive optics— instrumentation: spectrographs— planets and 
satellites: detection— techniques: high angular resolution 


1. INTRODUCTION 

Observing the population of planetary and brown 
dwarf companions orbiting young (~5-10 Myr) stars, 
soon after the dissipation of the primordial gaseous 
disk, is a key measure that will lend support to 
competin g formation models of substel l ar objects (e.g. 
Delgado-Donate et all 120041 IStamatellos fc Whitworthl 
20091) . Specifically, direct measurements of the orbital 
distribution of these objects shortly after formation (e.g. 
iDodson-Robinson et all 120091 : iKratter et all I2010D will 
serve as essential constraints to theoretical and numeri¬ 
cal models of planetary formation. Thus, o bserving low 
mass companions as early as possible (e.g. iKraus et all 
120141) will then serve as a “snapshot” of nascent system 
architecture, and largely eliminate any confusion about 
the initial conditions of companion for mation caused by 
subsequent dynamical proce sses (e.g. iScharf &; Menoul 
120091 IChatteriee et al.l [20Toh . Moreover, observing the 
luminosities of substellar companions in the first few 
million years is essential to constrain models of the in- 
tial entropy and temperatures of substellar objects (e.g. 
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iFortnev et al.ll2008t iMarlcau fe Cummindl2014ih 

However, due to the scarcity of young stars in the solar 
neighborhood, assembling a statistically robust sample 
of low-mass companions requires observations of large, 
newly-formed stellar associations. One such region, the 
Scorpius-Centaurus (hereafter “ Sco-Cen”) association 
with a distance of ~12 0-150pc (Ide Zeeuw et all 119991) 
and 5-20 Myr age (e.g. iPecaut, et al.l 12012b iSong et all 
I2012T) is the nearest OB Association. The young age of 
this association ensures that any planetary and brown 
dwarf companions will have elevated luminosity (e.g. 
iBaraffe et all 1200a) . allowing access to planetary mass 
objects even with observations achieving modest relative 
contrasts. 

In addition to selecting young targets, observing stars 
more massive than solar (“intermediate mass stars”, 2- 
5 M 0 ) may also enhance the probability of detection 
of substellar objects. Indeed, some studies suggest the 
fraction of super-Jupiter compa nions may be higher for 
stars more massive t han ~2 M^ (ICrenp fc J ohnsonll2011l: 
Vigan et al.; 2012 : iRameau et all 120131 : iReffcrt et all 
20151) . possibl y due to initially mor e massive circumstel- 
lar disks (e.g. lAndrews et aI1l2013l ). or possibly to serve 
as a reservoir for the conserved initial an gular momen¬ 
tum o f the star forming cloud (e.g. iKouwenhoven et al.l 
l2007b( ). The relatively short main-sequence lifetimes of 
these intermediate mass stars implies young associations 
such as Sco-Cen should have a greater fraction of these 
stars compared to the local solar neighborhood, making 
Sco-Cen a particularly promising region to study young, 
intermediate mass stars. 


To date, the only companions to Sco-Cen intermediate- 
mass stars lie at wide separa ti ons (~30-3 20 0 AU 


e.g. Kouwenhoven et al.l 2007bfc 

Aller et al.l 20131: 

Janson et al.l 2013E Lafreniere et all 

2014f Bailcv et all 

20141). The ~120-150pc distance 

to Sco-Cen stars 


means that orbital separations of substellar companions 
to these stars located near the ice line (5-10 AU), where 
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TABLE 1 

Table of Observations 


Target 

(HIP) 

Region 

^membership 

(%) 

SpT 

dist 

(pc) 

V 

(mag) 

Wl 

(mag) 

Ay 

(mag) 

■Mboi 

(mag) 

log(T eff ) 

(K) 

M 

(M 0 ) 

Observatory 
& UT Date 

71724 

UCL 

94% 

B8.5 

157±S 

6.63 

6.82±0.05 

0.04i0.23 

0.03i0.55 

4.072i0.087 

Q 41+0-53 

o.75 

VLT: 2011 June 10 

73990 

UCL 

92% 

A9 

125±i| 

8.23 

7.28±0.03 

0.30F0.17 

2.48i0.28 

3.872i0.022 

-i 7 Q+O.O 8 

VLT: 2011 June 10 

74865 

UCL 

96% 

F4 

115+i® 

9.00 

7.78±0.03 

0.20F0.07 

3.50i0.31 

3.822i0.006 

i 42+0-04 

VLT: 2011 June 10 

78196 

USco 

95% 

A0 

127±| 

7.03 

7.09±0.05 

-0.05i0.25 

1.36i0.45 

3.987i0.076 

9 4fi + 0,31 

^•^O_ 0 .59 

Keck: 2010 May 24 

78233 

USco 

79% 

F 0 

145+is 

—15 

9.17 

7.64±0.03 

0.84i0.11 

2.55i0.25 

3.857i0.012 

1 ^7+0.07 
i - D/ - 0.12 

Keck: 2010 April 25 

79124 

USco 

93% 

A0 

12311 1 

7.83 

6.96±0.05 

0.88i0.25 

1.30i0.46 

3.987i0.076 

9 ,4Q+0-32 
^• 4 °-0.58 

Keck: 2010 April 05 


Note. — “UCL” and “USco” refer to the Upper Centaurus Lupus, and Upper Scorpius regions of Sco-Cen, respectively. Membership probabilies 
in column three are taken from fRizzuto et al.l d201lP 


planet formation i s thought to be most efficient (e.g. 
iPollack et al.lll996f l. corresponds to angular separations 
very close to the near-infrared diffraction limit of 10 m 
telescopes ( X/D « 30-45 miliarcsec). Thus, sensitivity 
near, and within, the diffraction limit of large telescopes 
is needed to access smaller orbital separations for 
Sco-Cen stars (IKraus et alll2008l 1201 111 . . _ 

Aperture maski ng interferometry (e.g. iTuthill et al.1 
120001 : [Ireland 2013L and references therein), provides sen¬ 
sitivity at scales up to, and somewhat within, the usually 
defined diffraction limit ( /D — AX/D ~ 20-300 mas 
for Keck L' -band imaging). Appl i cations of this tech¬ 
nique (e.g.. Ilreland &; Kraud 120081 : IHinklev et al.l [2QTT1 : 
IKraus fc Irelandl I2012L and~references therein) use AO 
along with an opaque mask containing several holes, 
constructed such that the baseline between any two 
holes samples a unique spatial frequency in the pupil 
plane. Further, no coronagraphic mask is used, which 
avoids problems associated with measuring the relative 
astrometry between th e occulted h ost st ar and a de¬ 
tected companion (e.g. iDigbv et al.ll2006D . Despite its 
very good sensitivity to small inner working angles, aper¬ 
ture masking interfero metry usually only achieves typical 
contr asts of ~5-8 mag (IKraus et al.l 120081 IHinklev et al.l 
M). Nonetheless, these modest contrasts are still 
sensitive to planetary mass companions at y oung ages 
(Ilreland fe Krausll2008l : IKraus &; Irelandl[2012D . 

This paper presents discoveries of companions with ex¬ 
treme mass ratios, <?=M compan i on /Mi 1 ost Ri 0.01-0.08, from 
an ongoing multiplicity survey of ^140 intermediate mass 
stars (SpT=B0-F2) in the Sco-Cen region using aper¬ 
ture masking interferometry. We have selected our tar¬ 
gets based on the refined Bayesian Sco-Cen membersh ip 
selection technique described in iRizzuto et al.l (1201111 . 
which uses radial velocity information to confirm or re¬ 
ject candidates. Since high-mass sta rs frequent ly host 
one o r more binary companions (e.g. iDuchene fe Kraua 
I2013T) which prevents the required contrast from being 
achieved, the targets in this study have been screened us¬ 
ing past literature to eliminate binary systems with typi¬ 
cal separations of >30 mas and masses >0.1 Mg. Follow¬ 
ing this, we used results from our own ongoing snapshot 
AO imaging programs at Keck, VLT, and Palomar Ob¬ 
servatories to eliminate other systems with evidence for 
binarity. A more comprehensive discussion of our target 
selection and broad survey results will be given in a sub¬ 
sequent work (Rizzuto et ah, in prep). In ^ we describe 


the host star properties for our discoveries, followed by 
our observations and analysis strategy (10). In JHand (0 
we summarize our findings and place them in context. 

2. TARGET STAR PROPERTIES 

In TableQ] we list the basic properties of the targets 
described in this work. With an ultimate goal of calcu¬ 
lating host star mass, we start by calculating the bolo- 
metric magnitude (Mboi) from the V magnitude plus 
an estimate of the visual extinction Ay and bolometric 
correction BCy- We first estimated the visual extinc¬ 
tion for each of our targets by comparing the observed 
(V — K ) colors of our target stars with the (V — K ) 
colors for the corres ponding spectral types tabulated by 
iPecaut et al.l (1201211 . As noted in TableQ] the sp ectral 
types we assu me come either from the HD catalog dHoukl 
ri978h or from IPecaut fc Mamaiekl (I2013H . Our spectral 
types have uncertainties of ±2 subclass for the spec¬ 
tral typ es listed i n th e HD catalog, and ±1 for those 
listed in IPecaut et al.l ( 2012 ). We assume that our un¬ 
certainties in spectral type result in (V — K) color un¬ 
certainties of 0.1-0.2 mag, while the uncertainties in 
the color versus spectral type relation do not exceed 
0.05 mag. Next, using this estimated and the observed 
(V — K) color, we esti mate Ay- using the relation be¬ 
tween Ay and Ak from lSchlegel et aIT ( 1998 0. Using the 
spectral t ypes and uncer t aintie s, we adopt BCy values 
from IPecaut fc Mamaiekl (I2013H . where the uncertainty 
in BCy is set by the uncertainty in the spectral type. We 
then combine V mag, Ay, BCy, and the distance mod¬ 
ulus to calculate bolometric absolute magnitudes Mboi, 
and its uncertainty. Next, we estimate T e g and its uncer- 
tainty using the Sp T along with the tabulated values in 
IPecaut fc Mamaiekl (120131) . Lastly, we use the T e s and 
Mboi values calculated above to calculate the mass for 
each target star by calculating a two-dimensional sur¬ 
face, mass(T e g, Mboi), and then marginalizing the prob¬ 
ability distribution function of (T e g, Mboi) to obtain ler 
confidence intervals. These masses are listed in TableQ] 

As TableQ] shows, all of the stars considered in this 
paper are members of either the Upper Centaurus Lu¬ 
pus (hereafter “UCL”) or the Upper Scorpius (hereafter 
“U Sco”) subgroups. Rather than assign distances to 
each star that reflect the average distances of the Sco- 
Cen subgroups (e.g. 125±15 pc and 145±15 pc for UCL 
and U Sco, respectively), we use the individual Hippar- 
cos parallaxes recorded for these stars, since the indi¬ 
vidual parallax uncertainties associated with each target 
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TABLE 2 

Properties of Newly Detected Companions 


Companion 

A L' 
(mag) 

Sep 

(mas) 

PA 

(degrees) 

Confidence 

Level 

Mu 

(mag) 

10 Myr Mass 

(Mjup) 

20 Myr Mass 

(Mjup) 

Sep 

(AU) 

Q —A^comp/ Af s tar 
(10 Myr, 20 Myr) 

HIP71724B 

4 86+ 0 09 
‘i.oU —0 08 

10311 

20211 

>99.9% 

5 70+°- 33 

u 1 -0.34 

1461f 2 

21712 

1618 

0.0418-“, o.o6+8;81 

HIP73990B 

6 04+°- 52 

161±53 

28111 

>99.9% 

7 84+0.77 

4 -0.62 

21130 

501H 

2011 

°.°il8:S!. 0 03 lg.gl 

HIP73990C 

5 95 + 0 - 55 

^•^-0.44 

254+“ 

4811 3 

>99.9% 

7 75 +o.so 

4 • 4 -0.72 

22!f 

54118 

3211 

o.oil8;8?, 0.0318:88 

HIP74865B 

5.03«;f 3 

201 ±" 

11511 

>99.9% 

7 +0.59 

4 * O1 -0.60 

281?S 

65111 

2311 

o.o2i8;8i, o.o4!8;8i 

HIP78196B 

4.6l!g;i| 

74+10 

4 -10 

26611 

<99% 

6-1818:2 

98118 

152111 

918 

o.o4i8;8?, 0.0618:81 

HIP78233B 

4-72+g;“ 

13313 

201 j 

99.9% 

6-5518;“ 

861- 

1241H 

1918 

0.0518:81,0.0718.81 

HIP79124B 

4-30i°!° 

17711 

2421) 

99.9% 

5.8118:2 

1351H 

2011H 

2211 

0.0518:88, O.OS18.81 


Note. — Masses are derived from [Chabrier et al.l (120001) and IBaraffe et al.l (119981 . DoO.'.ji models. 


are comparable to the subgroup distance uncertainties 
(-10%). However, HIP 78233 has a recorded Hipparcos 
parallax inconsistent with the median U Sco members, 
with large uncertainty (4.84±1.37mas). So we assign to 
it a distance of 145±15 pc, the median distance to U Sco. 

While we have selected our targets partly based on 
their high proba bil ities of Sco-Cen membership as stated 
in iRizzuto et al.l (l2011f) . placement of the host stars on an 
HR diagram using the values calculated in Tabled] veri¬ 
fies all targets ca n be we ll fit be t ween the 10 and 20 Myr 
isochrones from iBressan et alJ (120121) . These ages are 
broadly consistent with the reported ages of UCL and U 
Sco. Some disagreement p ersists over t he age of U Sco , 
which ranges from 5 Myr (]Preibisch fc Zinnecker| |1999|: 
Preibisch et al.ll2002t ) to as high as 11 Mvr dPecaut et ail 
20121) . Indeed, a 5 Myr age is required to place low mass 
U Sco stars on an HR diagram (Rizzuto et a l. in prep , 
iPreibisch Sz Mamaieldl2008l : lKraus et al.ll2015f) . Further, 
when placed on an HR diagram at least half of our tar¬ 
gets (HIP71724, HIP78196 and HIP79124) have posi¬ 
tions in the HR diagram that are consistent with a 5 Myr 
isochrone. Nonetheless, for overall consistency we report 
our derived masses using a 10-20 Myr age range. 

3. OBSERVATION STRATEGY & ANALYSIS 

All the data presented in this work were obtained at 
L'-band wavelengths (3.76 ^irn) using the NIRC2 infrared 
camera and AO system at the W. M. Keck Observatory 
using , as well as the ESO Very Large Telescope (“VLT”) 
NACO AO system and infrared camera. An observing 
sequence consisted of observing a target star in two op¬ 
posed quadrants of the infrared camera: NIRC2 at Keck, 
and CONICA at the VLT. A nine hole aperture mask is 
used at Keck, and a seven hole mask was used at VLT, 
producing inte rferog rams like that shown in Figure 1 of 
iHinklev et ab (|2011l) . At both detector positions, we typ¬ 
ically obtained 15 images (30 images total) with an effec¬ 
tive exposure time of 20s each. Usually two to four such 
30-image sequences of each target star were obtained. 

We did not explicitly observe calibrator stars for each 
of our target stars, as is common practice. Rather, we 
use all of the stars in a given observing night as mu¬ 
tual calibrators. Those with closure phase signals in¬ 
dicative of a companion are weighted lower in the list 
of calibrators than those without. In a single observing 
night, this method allows roughly twice the number of 
Sco-Cen targets to be observed. To use the closure phase 
quantity to search for companions, we follow the analysis 


outl ined inlKraus et- al.l (|2008T ) , llreland fc Fraud (|2008f i 
and IHinklev et al.l (1201 111 , briefly summarized here. The 
data are initially flatfielded, sky subtracted, aligned, and 
corrected for cosmic rays. The bispectrum, the complex 
triple product of visibilities defined by the three baselines 
formed from any three subapertures, is then calculated. 
The phase of this complex quantity i s the closure phase. 

As d iscussed in lKraus et al.l (|2008l ) and IHinklev et all 
(|2011h . the calibrated object closure-phase is found by 
subtracting a weighted average of the closure-phase for 
the calibrator stars. For the analysis in this paper, which 
is motivated by the search for point sources, the squared 
visibilities were not used as they were noisier than the 
closure-phases. Each target was calibrated against all 
other stars to search for any deviations from single 
point-like sources. The Root Mean Square (RMS) cal¬ 
ibrated closure-phase was found for each of these target- 
calibrator pairs, and all calibrations that resulted in an 
RMS closure phase more than 1.5 times the minimum for 
each target over all calibrators were assigned a weight of 
zero. In practice, this meant that each target data set 
was calibrated by an average of the two or three calibra¬ 
tor data sets obtained closest in time. The lack of perfect 
closure phase calibration is still the dominant source of 
closure-phase noise in this analysis. 

4. RESULTS 

In Tabled we present the key properties of our newly 
detected companions, including the relative //-band con¬ 
trast ratios, angular separations, position angles, abso¬ 
lute magnitudes for the detected companions, as well 
as confidence levels of each of the detected compan¬ 
ions to these Sco-Cen BAF stars. For masses less than 
0.1 Mq, we convert these absolute magnitudes into physi- 
cal masses for the DUSTY models dChabrier et al.ll2000D 
corresponding to ages of 10 and 20 Myr. For masses 
greate r than 0.1 Mq, we interpolate the iBaraffe et all 
(1998:) models, assuming a mixing length parameter 
U m i X =1.9L/p, with H p denoting the pressure scale height. 
Our L'-band photometry for the host stars, used to cal¬ 
culate the relative brightness of the companions, is ob¬ 
tained from the WISE W1 channel. 

Figure]!] shows our newly detected objects in a mass 
versus semi-major axis diagram, where we plot the lower 
(upper) limit age of 10 Myr (20 Myr) for each object us¬ 
ing the values tabuated in Tabled For context, Figure]]] 
also shows several other detections of objects orbiting 
Sco-Cen stars of spectral type B, A, and F using conven- 


























































4 


Hinkley et al. 



Fig. 1.— Companions to Sco-Cen stars of spectral type B, A, and F (>1.5 Mq) expressed in terms of Jupiter masses ( Mj up ) and projected 
orbital separation (AU). The blue and white circles denote our new detections of companions to Sco-Cen BAF stars using aperture masking 
interferometry (see Tablc[2]) 5 assuming a lower limit age of 10 Myr (blue circular points), and an upper limit age of 20 Myr (white circular 
points). The gray symbols denote previous identifications from the literature of companions to Sco-Cen BAF stars obtained through 
conventional AO i maging (triangles, Kouwenhoven et al. 2005), interferometry (squares, Rizzuto et al. 2013) and more recent imaging 
studies (pentagons, Janson ct al. 2013; Lafrcniere ct al. 2014: Bailcv et al. 2014]). For compariso n, the brown circula r symbols indicate the 
planetary mass com panions to the young A-stars HR 8799, (3 Pic, and HD 95086, respectively (Marois ct al. 2010; Lagrange et al. 2010; 
I Rameau et al.|[2013D 
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tional AO (Kouwenhoven et al.l 2005 

; IJanson et al.l 2013t 

lLafreniere et all 20141: 'Bailev et al. 

20141). Also in the 


figure, we also show recent interferometric detections of 
com panions with primaril y stellar masses (>1 Mq) taken 
from lRizzuto ct all ( 2013 1. In addition, the Figure shows 
the directly imaged exoplanets orbiting the closer, young 
^10-30 Myr)A-stars HR 8799, /3 Pic, and HD 95086b 
Marois et al.ll20ldt lLagrange et al .1 [20 1 Ot I Rameau et al.l 

I2013T) . Our objects occupy a similar orbital range (10 

to 30 AU) to these directly ima ged exop lane ts, bu t with 
larger masses. Some works (e.g JVigan et al.l 120121 ) have 
suggested the peak of the companion distribution lies in 
this orbital range. 

Three of our newly detected objects, HIP 73990B, 
HIP73990C, and HIP74865B, have mass ratios clearly 
below the Hydrogen burning limit (~72 Mj up ), even 
assuming the older 20 Myr age. These are objects 
at 10-30 AU that unambig uousl y occupy the so-called 
“brown dwarf de s ert” (e.g. iKouwenhoven et ahl l2007at 
iKraus et al.ll2(M l20H . an observationally determined 
dearth of brown dwarf objects traditionally categorized 
as having g<0.1. Any detection of objects in this mass 
range will be particularly important to inform the oreti¬ 
cal a nd numerical models of multiple systems (e.g. iBatd 
120091 and references therein). 

Given the very small angular separations at which 
aperture masking interferometry performs (<0.25 ,, ) ; the 
likelihood of contamination from background stars is neg¬ 
ligible. For our survey, the number of expected con¬ 
taminants has been estimated based on the local sur¬ 
face density of stars in the vicinity of the targets in 
our sample. This density was estimated using both the 
2MASS survey and a star count algorithm that combines 
a spatial and luminosity model for the Milky Way (for 
the thin/ thick disk, halo, bulge and present-day mass 
function dReid et all 12002 ). Nonetheless, to obtain the 
number of expected false detections in our broad survey 
of 140 stars , we e xtrapolate the false alarm rate from 
IKraus et al.l (120081) . which predicts 0.3 false detections 
for 60 Sco-Cen stars. With this rate in hand, for our 
survey of 140 stars, we would expect ~0.7 false detec¬ 
tions. Thus, the probability of identifying seven com¬ 
panions, as we have done in this work, would require a 
false a larm rate an or der-of-magnitude greater than that 
of IKraus et al.l (I2008D . Such a circumstance is exceed¬ 
ingly unlikely, since the current study also focusses on the 
Sco-Cen region. Thus, the conventional need to confirm 
common proper motion is much less urgent. Nonetheless, 
we will continue to monitor these targets with the goal 
of fully characterizing the orbital motion of the compan¬ 
ions, as well as establishing common proper motion with 
the host stars. 

Of the six stars targetted in this study, only HIP 73990 
has evidence for significant excess emission at 22 /im as 
measured by the NASA WISE mission ([3.6]-[22]/xm = 
1.62), suggesting the presence of a debris disk. At the 
same time, none of our six target stars have any sta¬ 
tistically significant excess emission at 4.6 /im or 12 /xm, 
which would suggest the presense of more optically thick, 
protoplanetary disks. HIP73990 is also our only target 
with more than one detected companion. We devoted 
several additional post-processing tests to ensure the 
dual detections were not spurious (e.g. optical “ghosts” 
in the field of view, a companion present around an¬ 


other mutual calibrator star, etc.). A thorough repro¬ 
cessing of the data while varying the calibration scheme 
revealed that both detections were present irrespective 
of the number of mutual calibrators stars used. Further¬ 
more, the NACO L27 camera and L'-band filter have well 
characterized optical ghosts present in the focal plane. 
However, the position of these artifacts are fixed, and 
are easily excluded by a mask. Nonetheless, as an addi¬ 
tional check, we performed a data reduction using only 
those files with the interferograms in identical places on 
the focal plane. The detection of the companions were 
robust against these changes as well. Lastly, we rule out 
the possibility that the closure phase signal could per¬ 
haps be caused by aliasin g from a more distant source, 
since iJanson et all (|2013l) do not mention any detection 
of faint companions in their observations of this star us¬ 
ing conventional AO imaging. 

5. CONCLUSIONS & SUMMARY 

In this work, we report the detection of seven compan¬ 
ions to intermediate-mass stars of spectral types B, A, 
and F in the Sco-Cen association. With assumed ages 
of 10-20 Myr, our newly detected objects are observed 
shortly after the epoch of formation, residing on orbital 
scales comparable to objects in our solar system (~10- 
30 AU). The young age of these stars allows detection 
of brown dwarf (and potentially planetary) mass objects 
even with modest achieved contrasts AL'faA-6 mags. We 
highlight our main findings in this work as follows: 

1) All of our newly detected companions are unam¬ 
biguously “brown dwarf desert” objects with mass ratios 
g~0.01-0.08. The derived masses of our newly detected 
companions suggest they are more massive analogs of the 
planetary mass companions HR8799bcde, dMarois et all 
120101 ). HIP 73990 is perhaps the strongest analog for HR 
8799, with multiple detected objects with masses ~20-50 
Mj up (g«l-3%). Along with HIP 74865B, these three 
objects are all clearly below the Hydrogen burning limit, 
even assuming an age of 20 Myr. 

2) The objects presented in this paper have much 
smaller orbital separations than the previously reported 
substellar companions to S co-Cen BAF stars detected 
through conventional AO (IKouw enhove n et al.l l2007bl ; 
IJanson et al .11201:1 lLafreniere et al.ll2014D . 

3) The infrared brightnesses presented here will fur¬ 
ther serve as luminosity measurements of very young 
objects against which evolutionary models can be com¬ 
pared, thereby constraining the i nitial entropy of formin g 
low mass companions (e.g JMarlcau fc Cummin5l2014D . 

4) This study also begins to fill an important gap in our 
knowledge of the multipli city o f intermediate m ass stars 
at young ages (e.g. iDelgado-Donatc et al.ll2004D. Spec if- 
ically, some studies suggest (e.g. IKouwenhoven et all 
12007b! ) that multiplicity is an essential outcome of the 
formation of intermediate mass stars, serving as a reser¬ 
voir for the conserved initial angular momentum. 

5) The new companions presented in this paper will be 
prime targets not only for follow-up spectroscopy (e.g. 
iHinklcv ct al.l [20T5h . but also to search for fainter com¬ 
panions at large separations using the latest generation of 
dedicated exoplanet imagers such as GPI and SPHERE. 
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Finally, as GAIA parallaxes are derived for these objects, 
follow-up high resolution host star spectroscopy will be 
highly beneficial to better determine the host star phys¬ 
ical properties, such as T e g, log (g), and metallicity. 
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